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Abstract 

Numerical calculations of the kinetic energy of various extensions of the one- 
dimensional Hubbard model including dimerization and repulsion between 
nearest neighbours are reported. Using the sum rule that relates the kinetic 
energy to the integral of the optical conductivity, one can determine which pa- 
rameters are consistent with the reduction of the infrared oscillator strength 
that has been observed in the Bechgaard salts. This leads to improved esti- 
mates of the correlation parameters for both the TMTSF and TMTTF series. 
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The problem of finding an accurate description of the electronic properties of quasi-one- 
dimensional organic conductors has a long history. Even for the one-dimensional properties, 
which can be observed at not too low temperatures, there is no consensus. One of the origins 
of the difficulty is that one can work within two different kinds of model. The first one is the 
Fermi gas model also known as g-ology |T],|2| . This is the appropriate framework to describe 
the low-energy properties which are not of the usual Fermi liquid type but of the Luttinger 
liquid type. Most of the low-energy properties have been extensively analyzed within this 
type of model, and some information on the size of the parameters is available The 
second kind of model is the Hubbard model and its extensions. They are described by the 
Hamiltonian 

ieven,<T iodd,cr i i 

The parameters of this model are: i) a hopping integral t\ for the short bonds; ii) a hopping 
integral (< t\) for the long bonds; iii) an on-site repulsion U; iv) a repulsion between 
nearest neighbours V. In the following, energies will be measured in units of ti, and the 
basic dimensionless parameters are t<ijt\ for the dimerization and U/t\ and V/t\ for the 
Coulomb interactions. These parameters are large energy scales, and such a description is 
the natural framework to make contact with quantum chemistry or to analyze high-energy 
properties of the materials. 

Recently, it has been shown that, even for the low-energy properties, this description 
in terms of high energy scales can be very useful ||. In the context of the Luttinger 
liquid theory, the central parameter that describes the low-energy physics, the exponent 
K p , can take any value from to +00. However, for the quarter-filled, extended Hubbard 
model, which is the appropriate description of (TMTSF^PFg as far as quantum chemistry is 
concerned, it was possible to show that K p cannot be smaller than 1/4. It is actually possible 
to calculate K p as a function of the correlation parameters with a reasonable accuracy using 
standard numerical techniques anc ^ an accura t e determination of these parameters 

will be an important step toward a good understanding of the low-energy properties of the 
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Bechgaard salts. 

How can one determine the correlation parameters? In systems like transition metal 
compounds, the combination of photoemission and Auger spectroscopy has proved to be the 
most powerful tool to answer that question ||, but this is hopeless in the case of the Bech- 
gaard salts because the photoemission spectrum is itself the subject of much controversy 
H|5|]. An analysis of the charge-transfer spectra of many charge-transfer salts by Mazum- 
dar and Dixit led them to the general conclusion that U is about 1.5 eV and V is 0.5 eV 
10| . This analysis is certainly relevant, but in view of the very numerous experimental 



results that are now available for the Bechgaard salts, an analysis more specific to these 
compounds, including in particular the effect of dimerization, should be performed. The 
interpretation of the temperature dependence of the susceptibility |TT|j3[] and of the resis- 



tivity |]T2J] has given some preliminary information, but it is not very precise: Quantitative 
calculations are already very difficult for the Hubbard model with only on-site repulsion, 
and accurate estimates for the model of Eq. (1) as a function of t 2 /^i 5 U/ti and V/t\ are 
not available. Finally, quantum chemistry calculations of the correlation parameters have 
also been performed for these systems , but the reliability of the calculation is again very 
controversial. 

There is however one set of experimental data that has not been systematically used, 
namely the reflectivity measurements. That these results contain information on the local 



correlations has already been explained in great details by Jacobsen ||14|| . The idea is the 
following. On one hand, one can determine the plasma frequency from a Drude fit of the 
reflectivity spectrum. On the other hand, one can extract the optical conductivity through 
a standard procedure. Comparing the integral of the real part of the conductivity with the 
plasma frequency yields an estimate of the reduction of the kinetic energy due to correlations. 
To be more precise, let us denote by T the kinetic energy operator defined by 

T=-ti J2 ( c La+V + h - c - *2 E ( c *W liff + h.c.) (2) 

zeven,cr iodd,a 

The kinetic energy is given by E^m =< T >, where the expectation value is calculated in the 



ground-state of the full Hamiltonian H . The kinetic energy per site 4m is then defined as 
lim^ +00 Ek in /L, where L is the number of sites. The plasma frequency provides an estimate 
of the kinetic energy £° in calculated in the ground-state of T, i.e. without correlations, the 
integral of the conductivity provides an estimate of £ki n , and the reduction of kinetic energy 
is defined as the ratio £^i Q /S^ in . The main difficulty is where to stop in performing the 
integral of the conductivity. This can be a serious problem because, for strongly correlated 
systems, spectral weight coming from the conduction band can be found at high energy - 
typically around U - and it is impossible to disentangle this weight from other contributions 
to the conductivity. In the case of the Bechgaard salts, this is not too serious because the 
systems are essentially quarter-filled. In that case, it was shown by Maldague [fjl|, and 



recently confirmed by Eskes and Oles ||17|| , that the sum rule is almost exhausted by the 
lower band, and the estimates of the kinetic energy obtained by integrating over the lower 
band only are accurate to within a few percents. 

Performing such an analysis for the Bechgaard salts, Jacobsen reached the conclusion that 
£ k in/£kin is about °- 85 for (TMTSF) 2 C10 4 and 0.73 for (TMTTF) 2 PF 6 . Various estimates 
of the hopping integrals for (TMTTF) 2 PF 6 and (TMTSF) 2 C104 have been proposed on 
the basis of experimental results and quantum chemistry calculations. While there is still 
some uncertainty concerning their absolute value, especially in the case of (TMTTF) 2 PFg, 
the ratio t 2 /ti is believed to be approximately equal to 0.9 for (TMTSF) 2 C10 4 and 0.7 



for (TMTTF) 2 PF 6 [ 15(1 ■ So to extract information about the correlation parameters from 
Jacobsen measurements, one just needs accurate estimates of the kinetic energy in the 
ground-state of the Hamiltonian of Eq. (1) as a function of U/ti and V/t\. When Jacobsen 
published his results, nothing of the sort was available, and he could not go beyond a 
qualitative analysis of the results based on the numerical study of a system of 2 particles 
on 4 sites. Motivated by Jacobsen's results, Baeriswyl et al [18|] calculated the kinetic 
energy of the standard Hubbard model using the Bethe ansatz solution at half-filling and 
the Gutzwiller ansatz away from half-filling. Their results confirm the trends, namely that 
correlations have to be invoked to explain the reduction of oscillator strength, but they do 



not allow for a precise interpretation of the experimental results. We are not aware of any 
further work on that problem. 

In this paper, we calculate the kinetic energy finite clusters on the basis of numerical 
results obtained on finite clusters by exact diagonalization. For a given cluster, the ex- 
pectation value of the kinetic energy is most easily obtained using the Hellmann-Feynman 
theorem as 

-Ekin = h d tl E G .s. + h d t2 E G .s. (3) 

where the ground-state energy Eq.s. =< H > is evaluated with Lanczos algorithm. £kin is 
then obtained from a finite-size scaling analysis of E kin / L. For non-interacting electrons, 
it is easy to show that the finite-size corrections go as 1/L 2 , where L is the number of 
sites. This remains true for the total energy per site for Luttinger liquids. So it is quite 
natural to assume that this is also true for the derivatives of this quantity with respect to 
the hopping integrals, and thus for the kinetic energy. Our numerical results clearly support 
this assumption. A typical example is given in Fig. (1), where we have plotted E kin /L as 
a function of 1/L 2 for L — 8, 12 and 16. The two curves have been obtained for different 
boundary conditions correponding to closed and open shells respectively. The scaling law 
is quite accurately satisfied, and the two boundary conditions give estimates in very good 
agreement, which lends further support in favour of the assumption that the scaling is in 
1/L 2 . In most cases where we have tried both types of boundary conditions, the slope was 
smaller for the boundary conditions corresponding to open shells, and all the results given in 
the rest of the paper have been obtained with such boundary conditions. Besides, we have 
compared the results obtained by using only L = 8 and L = 12 with estimates obtained 
using also the results for L = 16 for a few cases, and the error was always less than 1%. So, 
unless one needs a very accurate value of the kinetic energy, it is sufliccient to use systems 
with 8 and 12 sites to perform the 1/L 2 extrapolation. This has been systematically done 
in the following. 

Let us start with the results obtained in the case V/t± — (no repulsion between first 
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neighbours). The ratio S^in/S^is plotted as a function of U/ti for three values of £2/£i in 
Fig. (2) for a quarter-filled system. In all the cases, it decreases with U, in agreement with 
the intuitive idea that correlations make the motion more difficult and lead to a reduction of 
kinetic energy. What is maybe more surprising is that, even for very large values of U/ti, the 
reduction is not so big. That our results are still valid for large U/ti can actually be checked 
quite easily. The value for U/t\ — > +00 is the same as for spinless fermions at half-filling 
(see Table I), and our results plotted as a function of t\/U extrapolate nicely toward that 
Uniting value (see Fig. (3)). So, taking for granted that £2/^1 = -7 for (TMTTF^PFg , a 
reduction of 0.73 is incompatible with on-site repulsion only. This can be considered as a 
proof that repulsion on neighbouring sites is important in these systems. 

Let us now consider the general case described by Eq. (1). The question we would like to 
answer is the following: What are the values of U /t\ and V jt\ that are compatible with the 
known values of £ 2 /£i and of S^/S^ for (TMTSF) 2 C10 4 and (TMTTF) 2 PF 6 ? The most 
convenient thing to do is to plot the curves of constant ^kin/^km m the {U/t\, V/ti) plane 
for the values of £2/^1 of interest. Such plots for quarter- filled systems are given in Fig. (4) 
for £2/^1 — 1, 0.9 and 0.7. The basic features of these curves are again quite natural. The 
only one that deserves a special comment is the reentrant behaviour for large V/t\ and small 
U jt\. Another way of looking at the same effect is to notice that, for a given value of V/ti, 
the kinetic energy first increases before it decreases for U/ £1 large enough. This presumably 
comes from local pairs, which are known to exist and to be very heavy objects in the small 
U/ti, large V/t\ limit ||19|| , and which become lighter when U ' jt\ increases. 



There is now no problem to find parameters that give a reduction of kinetic energy 
of 0.85 for (TMTSF) 2 C10 4 and of 0.73 for (TMTTF) 2 PF 6 . The corresponding curves of 
possible pairs (U/ti, V/ti) are reproduced in Fig. (5). This is not very useful unless we can 
decide where the actual parameters are located on these curves. This is actually possible on 
the basis of general arguments coming from quantum chemistry. Unlike the (/-parameters 
of the Fermi-gas model, U and V have a simple microscopic meaning: U is the energy 
needed to put two particles on the same site, and V is the energy needed to put them 



on neighbouring sites. Now, the molecules TMTSF and TMTTF are very similar, the only 
difference being that the 3p-orbitals of sulfur in TMTTF are more concentrated than the 4p- 
orbitals of selenium in TMTSF, and quantum chemistry calculations predict that the ratio 
U (TMTTF) /U (TMTSF) is in the range 1.0 - 1.25. We also know from different sources 
0,|T§ that t\ (TMTTF) /t x (TMTSF) is in the range 0.7 - 0.8. So the ratio U/t x is at most 
80% larger for TMTTF than for TMTSF. But from Fig. (5), we know that this ratio is 
at most 8 for TMTSF, so it is at most 15 for TMTTF. Looking again at Fig. (5), this 
means that V/t\ is at least equal to 2 for TMTTF. Now, on the basis of quantum chemistry, 
V is expected to be roughly the same for both types of compounds, which implies that 
V/t% is at least 1.5 for TMTSF. But according to Fig. (5), this means that U jt\ cannot 
be larger than 6 in that compound. This again puts a constraint on U/t\ in TMTTF, and 
so on. Finally, if we use the constraints given by quantum chemistry, we end up with the 
parameters given in Table 2. These values are just estimates, and one should put error bars 
on them. The main source of uncertainty probably comes from the experimental results. In 
particular, the saturation value of the integral of the loss function generally gives numbers 
in good agreement with the plasma frequency deduced from a Drude fit of the reflectivity 
spectrum. However, in the present case, the values obtained from the loss function are 
slightly larger 0, leading to smaller values of £ kin /££i n (0.67 for (TMTTF) 2 PF 6 , 0.80 for 
(TMTSF) 2 C104). Another source of uncertainty lies in the values used for the ratios t^/ti, 
but the dependence of 4in/ '^kin on this parameter is smooth (see Fig. 4). Finally, the 
location on the curves is only approximative because the arguments derived from quantum 
chemistry are only qualitative or, at best, semi-quantitative. 

In spite of these sources of uncertainty, two points seem to be clearly established. First, 
the reduction of oscillator strength reported for TMTTF implies that the repulsion between 
first neighbours is not negligible. Second, U/ti is larger in TMTTF than in TMTSF. If that 
was not the case, V would have to be much larger in TMTTF than in TMTSF, which can 
be rejected on the basis of quantum chemistry. This last conclusion should be contrasted 
with the interpretation of the temperature dependence of the susceptibility of Wzietek et 
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al [[|, which lead them to conclude that U/t\ is about the same in both series. Recent 
calculations of the susceptibility |TT| suggest however that the temperature dependence 
is actually consistent with our present conclusion that U/t\ is larger in TMTTF than in 
TMTSF. 

In conclusion, the reduction of oscillator strength reported for the Bechgaard salts by 
Jacobsen leads to precise and useful information on the size of the correlation parameters 
in these compounds. The main difference with respect to the microscopic models used so 
far in the interpretation of various experimental results (susceptibility, minimum of resis- 
tivity,...) lies in the presence of a relatively large value of the repulsion between nearest 
neighbours. Whether the estimates proposed in the present paper are consistent with the 
other experimental data remains to be seen. 

I am very grateful to Jean-Paul Pouget, who encouraged me to look at the optical data 
of Jacobsen. 



S 



REFERENCES 

[1] J. Solyom, Adv. Phys. 28, 201 (1979). 

[2] V. Emery, in Highly Conducting One-Dimensional Solids, edited by J. T. Devrese et al 
(Plenum, New York, 1979), p. 327. 

[3] C. Bourbonnais, J. Phys. (Paris) I 3, 143 (1993). 

[4] P. Wzietek, F. Creuzet, C. Bourbonnais, D. Jerome, K. Bechgaard and P. Batail, J. 
Phys. (Paris) I 3, 171 (1993). 

[5] F. Mila and X. Zotos, Europhys. Lett. 24, 133 (1993). 

[6] H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990); Int. J. Mod. Phys. B 5, 57 (1991). 
[7] M. Ogata et al, Phys. Rev. Lett. 66, 2388 (1991). 

[8] G. A. Sawatzky in High Temperature Superconductivity, eds. D. P. Tunstall and W. 
Barford, published by Adam Hilger, Bristol (1991). 

[9] B. Dardel et al, Europhys. Lett. 24, 687 (1993). 

[10] S. Mazumdar and S. N. Dixit, Phys. Rev. B 34, 3683 (1986). 

[11] F. Mila and K. Penc Phys. Rev. B 51, xxxx (1995). 

[12] K. Penc and F. Mila, Phys. Rev. B 50, 11429 (1994). 

[13] A. Fritsch and L. Ducasse, J. Physique I 1, 855 (1991). 

[14] See C. S. Jacobsen, J. Phys. C 19, 5643 (1986), and references therein. 

[15] L. Ducasse et al, J. Phys. C 19, 3805 (1986). 

[16] P. F. Maldague, Phys. Rev. B 16, 2437 (1977). 

[17] H. Eskes and A. M. Oles, Phys. Rev. Lett. 73, 732 (1994). 

[18] D. Baeriswyl, J. Carmelo, A. Luther, Phys. Rev. B 33, 7247 (1986). 

9 



[19] K. Penc and F. Mila, Phys. Rev. B 49, 9670 (1994). 



10 



FIGURES 

FIG. 1. Finite-size scaling of the kinetic energy for t 2 /h = 1, U/t\ = 10, and V/t\ = 0. Upper 
curve: closed shell; lower curve: open shell. 

FIG. 2. Kinetic energy as a function of U/t\ for V/t\ = 0. a) £2/^1 = 1; b) t 2 /ti = 0.9; c) 

t 2 /h = 0.7 

FIG. 3. Large U behaviour of the kinetic energy for V/t\ = 0. a) t 2 /t\ = 0.9; b) £2/^1 = 0.7 

FIG. 4. Constant kinetic energy plots from £ k in/£kin = °- 99 (bottom left) to ^kin/^kin = °- 50 
(top right), a) t 2 /h = 1; b) t 2 /h = 0.9; c) t 2 /h = 0.7 

FIG. 5. Comparison of the curves giving the correct reduction of kinetic energy. Crosses: 
(TMTSF) 2 C10 4 (t 2 /h = 0.9, S kin /£^ n = 0.85); Squares: (TMTTF) 2 PF 6 (t 2 /h = 0.7, 

^kin/4°in = 0.73) 
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TABLES 



t 2 /h 


1. 


0.9 


0.8 


0.7 


0.6 


0.5 


0.4 


0.3 


0.2 


0.1 


0. 


^kin/^kin 


.707 


.711 


.720 


.733 


.752 


.777 


.808 


.845 


.889 


.940 


1. 



TABLE I. Reduction of kinetic energy for U/ti — ► +oo and V/t\ = 



Compound t 2 /^i U/t\ V/t\ 



(TMTSF) 2 C10 4 


0.9 


5.0 


2.0 


(TMTTF) 2 PF 6 


0.7 


7.0 


2.8 



TABLE II. 
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